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The MAGIC collaboration reports the detection of the blazar S5 0716+714 
(z = 0.31 ± 0.08) in very high energy gamma-rays. The observations were per- 
formed in November 2007 and in April 2008, and were triggered by the KVA 
telescope due to the high optical state of the object. An overall significance 
of the signal accounts to S — 5.8 a for 13.1 hours of data. Most of the signal 
(S = 6.9 a) comes from the April 2008 data sample during a higher optical state 
of the object suggesting a possible correlation between the VHE 7-ray and optical 
emissions. The differential energy spectrum of the 2008 data sample follows a 
power law with a photon index of T — 3.45 ± 0.54 s ^ a ^ ± 0.2 S y S ^, and the integral 
flux above 400 GeV is at the level of (7.5 ± 2.2 stat ± 2.3 syst ) x 10" 12 cm" 2 s _1 , 
corresponding to a 9% Crab Nebula flux. Modeling of the broad band spectral 
energy distribution indicates that a structured jet model appears to be more 
promising in describing the available data than a simple one zone synchrotron 
self-Compton model. 

Subject headings: gamma rays: observations — BL Lacertae objects: individual 
(S5 0716+714) 



Introduction 



Blazars, a common term used for flat spectrum radio quasars (FSRQ) and BL Lacertae 
objects, appear to be the most energetic types of Active Galactic Nuclei (AGN). In these 
objects the dominant radiation component originates in a relativistic jet pointed nearly 
towards the observer. The double-peaked spectral energy distribution (SED) of blazars is 
attributed to a population of relativistic electrons spiraling in the magnetic field of the jet. 
The low energy peak is due to synchrotron emission and the second, high energy peak is often 
attributed to inyerse Compton scattering of low energy photons in leptonic e mission models 



(jMaraschi et al.lll992l ; iDermer fc Schlickeiser Ill993l : iBloom fc Marscherlll996l ). Models based 



Mannheim 


1993; 


Miicke et al. 


2003) 



For most FSRQs and a large fraction of BL 



submillimeter and optical. On the other hand, for most of the sources detecte d to emit 
VHE 7-rays (HBLaj) the low energy peak is located at UV to X-rays energies ( Padovani 



20071 ). The high energy peak is typically at MeV-GeV energies. Blazars are highly variable 



X LBL = low frequency peaking BL Lacertae 
2 HBL = high frequency peaking BL Lac 
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in all wavebands and the relation between variability in different bands is a key element in 
discriminating between different models. 

The MAGIC Collaboration is performing Target of Opportunity observations of sources 
in a high flux state in the optical and/or X-ray band. Optically triggere d observations have 
resulted in the discov ery of VHE 7-rays from Mrk 180 jAlbert et alJbood ) and 1ES 1011+496 
( lAlbert et al.ll2007al ). In this paper we report the results of observations of S5 0716+714 in 
November 2007 and April 2008. The observation at the latter date re sulted in the discovery 
of VHE 7-rays from the source as announced in iTeshima et al.l (120081 ) . 



The BL Lac object S5 0716+714 has been studied intensively at all frequency bands. It 
is highly variable with rapid variations observed from the radio to X-ray bands (Wagner et al. 
1996). It has therefore been target to several multiwavelength campaigns, th e most recent 



one organized by the GLAS T- AGILE Support Program in July-November 2007 (jVillata et al. 



20081 ; iGiommi et al.ll2008al ). Due to the very bright nucleus, which outshines the host galaxy, 
the reds hift of S5 0716+714 is still uncertain. The recent photometric detection of the host 
galaxy jNilsson et al.ll2008| ) suggests a redshift of z = 0.31 ± 0.08 which is consistent with 
the redshift z = 0.26 determined b y spectroscopy for three galaxies close to the location 
of S5 0716+714 Jstickel et al.lll993h . In the SED of S5 0716+714 the synchrotron peak is 
located in t he optical band and i s, therefore, classified either as LBL (INieppola et al.ll2006l ) 
or as IBlfl f lPadovani et~al"1ll995h 



S5 0716+714 was detected in the MeV energy range several times at diff erent flux levels 



by th e EGRET detector on board the Compton Gamma-ray Observatory (IHartman et al. 
19991 ). In 2008 AGILE reported the detection of a variable 7-ray flux with a peak flux 



density above the maximum reporte d from EGRET (|Chen et al.ll2008l ). S5 0716+714 is also 
on the Fermi-LAT bright source list (lAbdo et al.ll2009l ). Observations at VHE 7-ray energies 
by HEGRA resulted in an upper limit of F(> 1.6 TeV) = 3.13 x 10 -12 photons/cm 2 /s 
(lAharonian et al.l 120041 ) . In this paper we present the first detection of VHE 7-rays from 
S5 0716+714. It is the third optically triggered discovery of a VHE 7-ray emitting blazar 
by MAGIC. 



2. Observations 

The MAGIC (Major Atmospheric Gamma-ray Cherenkov) telescope is a standalone 
imaging atmospheric Cherenkov telescope located on the Canary Island of La Palma. MAGIC 



3 IBL = intermediate frequency peaking BL Lacertae 
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has a standard trigger threshold of 60 GeV for observations at low zenith angles, an angular 
resolution o f ~ 0.1° on the even t by event basis and an energy resolution above 150 GeV of 
-25% 



sec 



Albert et al.ll2008aL for details) 



The Tuorla blazar monitoring program^ ( Takalo et al. 2007 ) monitors S5 0716+714 
on a nightly basis using the KVA 35 cm telescope at La Palma|j and the Tuorla 1 meter 
telescope. At the end of October 2007 (22th, MJD=54395) the optical flux had more than 
doubled (from 19 mJy to 42 mJy) in less than a month and, according to a predefined criteria 
MAGIC was alerted. Due to moon and weather constraints, the MAGIC observations started 
11 days later, when the optical flux had already decreased significa ntly (see Fig. [Q. The 
observations were performed in the wobble mode (IDaum et al.lll997l ) pointing 0.4 deg offset 
from the source to enable simultaneous estimation of the background from the same field of 
view. MAGIC observed the source during 14 nights for a total of 16.8 hours. During some 
nights the observing conditions were rather poor and the affected data were rejected from 
the analysis. The exposure time for good quality data amounts to 10.3 hours. The zenith 
angle range of these observations was from 42 to 46 degrees. 

In April 2008 a new bright and fast optical flare occurred. The optical flux almost 
doubled within three nights (14th of April, MJD=54570: 29 mJy, 17th April: 52 mJy), and 
at 17th of April reached its historical maximum value. MAGIC started the observations 5 
nights later, when the moon conditions allowed. The source was observed during 9 nights 
with zenith angles from 47 to 55 degrees for a total of 7.1 hours. Unfortunately, during 
the last 6 nights of the observations there was strong calima wind carrying fine sand from 
Sahara desert and these data were, therefore, of bad quality and rejected from the analysis. 
The total exposure time of good quality data for this observation period amounts to only 
2.8 hours. The total live time of S5 0716+714 MAGIC observations in 2007 and 2008 after 
data quality cuts was 13.1 hours. 



3. Data Analysis and Results 



The MAGIC da ta were analy z ed us ing the standard analysis chain as described in 
Albert et al.l f|2008al lbh : | Aliu et al.l (120091 ). In order to suppress the background showers 
produced by c harged cosmic rays, a multivariate classification method known as Random 
Forest is used ( lAlbert et al.l l2008d ) . For every event, the algorithm takes as input a set of 



4 http://users.utu.fi/kani/lm 
5 http://tur3. tur.iac.es/ 
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Fig. 1. — The light curve of S5 0716+714 as measured from November 2007 until April 2008. 
The day-by-day 7-ray light curve from MAGIC is shown in the upper panel for the 2007 
data (left panel) and 2008 data (right panel), whereas the optical KVA data are shown in 
the lower panel. The simultaneous optical with the MAGIC data are marked black. The 
error bars (ler) of the optical fluxes are smaller than the points and thus not visible. 



image parameters, and produces one single parameter as output, called Hadronness. The 
background rejection is then achieved by a cut in Hadronness, which was optimized using 
Crab Nebula data taken under comparable conditions. 

The cut in | Alpha | that defines the signal region was also optimized in the same way. An 
additional cut removed the events with a total charge of less than 200 photoelectrons (phe) 
in order to assure a better background rejection. We used standard cuts in Hadronness 
and |Alpha|, which are determined to give the best significance for a point-like source with 
a flux on the 10% level of the Crab Nebula flux. For the given cuts and the relatively 
large zenith angles of the observations the analysis threshold corresponds to 400 GeV. The 
resulting |Alpha| distribution after all cuts for the overall S5 0716+714 data sample in 
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Fig. 2. — |Alpha| distribution after all cuts for the total MAGIC data sample in 2007 - 
2008. A 7-ray excess with a significance of 5.8 a is found. 

2007 - 2008 is shown in Fig. [2J An overall excess of 25 2 7-like events corresponding to a 



significance of S = 5.8 a was found (following Eq. 17 in [Li fc Mai (11983 ). N on = 1548, N ofT 
= 3887, alpha = 0.33). Most of the signal comes from the 2008 data sample: the analysis of 
the 2008 data alone results in 176 excess events over 422 background events corresponding 
to S = 6.9 a. From the 2007 data alone an excess corresponding to S = 2.2 a was found. 

The day- by-day light curve as measured by MAGIC data is shown in Fig. [1] (upper panel) 
together with the optical KVA light curve (lower panel). In November 2007 the MAGIC flux 
above 400 GeV is at F 2 oo 7 (> 0.4TeV) = (0.8 ± 0.7 stat ± 0.2 syst ) x lO-^cnrV 1 ], whereas 
the flux is about 9 times higher in 2008: F 200 8(> 0.4TeV) = (7.5 ± 2.2 stat ± 2.3 syst ) x 
10 [cm~ 2 s -1 ]. No significant variability is seen on time scales shorter than 6 months. 
Given the limited effective exposure times these observations are not sensitive to variability 
on shorter time scales. The individual MAGIC points would have low significance, and an 
intra-night variability by a factor of at least ten would have been required to detect it. In 
the optical band, instead, a clear variability on time scales from days to months is visible 
with two distinct flares: the first in October 2007, and the second in April 2008 (Fig. [1] 
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lower panel). 



The differential energy spectrum is c alculated only for th e April 2008 data set. The 
measured and unfolded for detector effects ([Albert et al.ll2007d ) spectrum is shown in Fig. [3J 
The data point at E — 1.3 TeV has a significance below la and was, therefore, converted 
into an upper limit corresponding to a 95% confidence level. The measured spectrum can 
be well fitted by a simple power law (with the differential flux given in units of TeV -1 cm -2 



diV 



(2.4 ±0.8) x 10~ n (£/500GeV) 



-3.5±0.5 



dE dA dt 

The errors are statistical only. The systematic uncertainties are estimated to be 0.2 on the 
photon index and 30% on the absolute flux level. Due to the energy- dependent attenua- 
tion of VHE 7-rays wit h low-energy photons of the extragalactic background light (EBL, 
Gould fc Schrederlll967l ). the VHE 7-ray flux of distant sources is significantly suppressed. 
W e calculated the de a bsorb ed, i.e. intrinsic, spectrum of S5 0716+714 using an EBL model 
of iFranceschini et al.l (120081 1 and assuming a redshift of z — 0.26. The resulting intrinsic 
spectrum (shown in Fig. red points) has a fitted photon index of Y = 1.8 ± 0.6, which is 
well within the range of other extragalactic sources measured so far. 

As the source redshift is still uncertain, we used the MAGIC spectra to calculate upper 
limits to the redshift. We assumed two different maximum values for a possible hardness of 
the intrinsic spectrum: 1.5, being a canonical value for a 7-ray spectrum emitted by electrons 
with a spectral index of 2.0; and 0.666, being the limiting case for a 7-ray s pectrum emitted 
by a m onoenergetic electron distribution. Using the method described in iMazin fc Goebel 
( 120071 ) we obtain the following upper limits for the redshift: z < 0.46 (for an assumed 
intrinsic spectrum with a power law photon index of 1.5) and z < 0.59 (for an assumed 
intrinsic spectrum with a power law photon index of 2/3). Both limits agree with the redshift 
determined from the host galaxy detection (z = 0.31 ± 0.08) and from the spectroscopy of 3 
nearby galaxies (z = 0.26). 



4. Discussion 

MAGIC observed the blazar S5 0716+714 in November 2007 and April 2008, the obser- 
vations resulting in the discovery of a very high energy 7-ray excess with a significance of 5.8cr. 
During the November 2007 MAGIC observations the average optical flux was ~ 20mJy, while 
in April 2008 the optical flux was ~ 45mJy. The same trend is also visible in the MAGIC 
data: the flux in April 2008 is significantly higher than in November 2007. This seems 
to support the indication seen in previous MAGIC observations for other BL Lac objects 
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Albert et al.l f|2006l . l2007al jbh . that there is a connection between optical high states and VHE 
7-ray high states. 



In April 2008 S5 0716+714 was also in a historical high state in X-rays (jGiommi et al. 



2008bj) and the optical polar ization angle started to rotate immediately after the optical 
maximum had been reached (jLarionov et al.l 120081 ) . However, the radio flux at 37 GHz did 
remain at a quiescent level (A. Lahteenmaki, priv.comm). This energy dependent behavior 
is very similar to the one seen in BL Lac in 2005, whic h was attributed to an emission feature 
moving in a helical path upstream of the VLB A core ( IMarscher et al.ll2008l ). The high VHE 
7-ray flux of S5 0716+714 observed by MAGIC might originate from such a moving emission 
feature, but a more detailed study of the multiwavelength light curves from spring 2008 are 
needed to confirm this. 

We tried modeling the SED with a one-zone synchrotron self Compton (SSC) model 
with relativistic electr ons following a smoothed broken power law energy distribution (see 
Tavecchio et al.l (120011 ) for a full description). The fit parameters are: the minimum electron 
energy 7 m i n = 10 4 , the break energy in the electron spectrum 7^ = 1.9 • 10 4 , the maximum 
electrom energy 7 max = 7 • 10 5 , the electron spectral index below and above the break 
energy n\ = 2, n 2 = 4.5, respectively, the magnetic field B = 0.13 G, the normalization 
factor K = 1.4 ■ 10 5 particles/cm 3 , the radius of the emitting region R = 7 ■ 10 15 cm and 
the Doppler factor 5 = 30. To account for the observed TeV flux corrected for intergalactic 
absorption the model (continuous line in Fig. H]) predicts a high intensity peak in the 10 GeV 
range. Although the source was flaring in optical at the time of the MAGIC detection, 
such high 7-ray flux (more than 10 times larger than observed by EGRET and Fermi, open 
blue bow-tie in Fig. H]) appears somewhat implausible. Therefore, we also considered the 
emission from a structured je t, modeled as a fast "spine" surrounded by a slower moving 
"layer" (iGhisellini et al.ll2005l ). Photons emitted by the layer are subject to IC scattering 
by the relativistic electrons of the spine thus introducing an additional contribution to the 
IC spectrum at energies higher than SSC can reach. The computed model is shown as a 
dashed line in Fig. HI The GeV excess is substantially reduced and the spectral shape in 
the 7-ray domain is consistent with the spectral indices measured at both GeV and TeV 
energies, though with presumably different intensity states. 

The discovery of high flux of VHE 7-rays by MAGIC from the blazar S5 0716+714 
in April 2008 should be compared to detailed light curves from other wavelengths (radio, 
optical and X-rays) as well as VLBA maps and optical polarization light curve in order to 
further investigate the origin of this interesting event. 
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Fig. 3. — The differential energy spectrum of S5 0716+714. Only data from April 2008 
are used. The measured (black points) as well as deabsorbed (using the EBL model of 



Franceschini et al.l ((20081) and assuming z=0.26, red points) spectra are shown. The results 



from a power law fit to the measured spectrum are shown in the plot. The Crab Nebula 
spectrum is shown for comparison (green dashed line). 
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Fig. 4.— The broad band SED of S5 0716+714. In red is the KVA (from 23 april 2008), 
Swift/UVOT and Swift /XRT (from 29 April 2008) and the deabsorped MAGIC (from 22 to 
24 April 2008) data. The MAGIC measured flux for April 2008 is shown in black and the 
MAGIC upperlimi t for November 2007 is shown with grey arrow. The grey symbols show 
historical data (see lTavecchio fc Ghisellini 1120091 ) . The bow-tie shows the flux from the Fermi 
bright list. The solid line shows the overall emission calculated with a one-zone SSC model 
using the April 2008 data (red), see text for model paramet ers. The dashed line sh o ws th e 
emission calculated with the spine- layer model described in iTavecchio fc Ghisellini I (120091 ). 
We assume the same parameters reported there, but a volume three times larger. 



